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Abstract. We use the Hubble Ultra Deep Field (UDF) to study the galaxy 
luminosity-size (M-R e ) distribution. With a careful analysis of selection effects 
due to both detection completeness and measurement reliability we identify bias- 
free regions in the M-R e plane for a series of volume-limited samples. We also 
investigate the colour-log(n) distribution of these galaxies and further subdivide 
our data by structural type to separately study compact and diffuse objects. By 
comparison to the nearby Millennium Galaxy Catalogue, we present tenative 
evidence for evolution of diffuse, disk-like galaxies with redshift — both in mean 
surface brightness and the slope of the M-R e relation. In contrast we find no 
evidence of structural evolution in the compact galaxy M-R e relation over this 
redshift range, although there is a suggestion of colour evolution. We also high- 
light the importance of considering surface brightness dependent measurement 
biases in addition to incompleteness biases. In particular, the increasing, sys- 
tematic under-estimation of Kron fluxes towards low surface brightnesses may 
cause diffuse, yet luminous, systems to be mistaken for faint, compact objects. 

1. Introduction 

Traditionally studies of galaxy evolution for the entire galaxy population have 
been pursued by looking for variations in the m onovariate lumin osity func- 
tion, either for the whole population (e.g., CFRS; [Lilly et aflll995), or subdi- 



Ilbert et al.l 



vided by morphological type or colour cuts (VVDS/COMBO-17; 
2006). However the resolution attainable with the Hubble Space Telescope 
enables one to move bey ond this to study multivariate distr ibutions, such as 



enables one to move bey ond tnis to study multivariate distr ibutions, sucn a; 
the luminosity-s ize plane (Shen et al.ll2003l ; [Driver et aill2005l ). and the colour 



structure plane ( Driver et al. 20061 ). These provide significantly more stringent 



constraints and are also now bec oming the focus of numerical simulations (e.g., 
lAlmeida. Baugh and Lacey| [2~006V In this article we summarise work under way 
([Cameron Driverl l2007bl ) to investigate both of these multivariate distribu- 
tions through the comparison of galaxies drawn from the Ultra Deep Field to 
our zero redshift reference sample defined by the MGC. We pay particular at- 
tention to the multiple selection biases inherent in this work. 
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Data 



We use the Hubble Ultra Deep Field (s eelBeckwith et al. 
ies out to z = 1.5 and the MG C (see iLiske et al.]2003l: 



2006 ) to sample g alax- 



Driver et al 



2005) for 



a z 



reference sample. In ICameron & Driverl (|2007al ) we present iUDF 
BRIGHT, a sample of 2497 sources brighter than 28th magnitude in the ACS 
UDF i-band image. These were detected by running SExtractor with a constant 
surface brightness threshold of 27.395 mag arcsec -2 . Kron magnitudes and ellip- 
tical aperture half light radii were computed for all these objects. Photometric 
redshift estimates based on B through //-band aperture-matched magnitudes 
and global Sersic ind ices for 2172 ga laxies in iUDF-BRIGHT were obtained 



from the catalogues of lCoe et al 



(2006). These redshifts and broadband magni- 
tudes were used to derive individual galaxy K-corrections (observed ACS i-band 
to rest-frame M GC 5-band) via fitting of SED templates from the library of 
PoggiantH ( 19971 ). We also included a flat spectrum model to approximate the 
K-corrections relevant to the extreme blue galaxies found in this type of deep 
imaging. The iUDF-BRIGHT K-corrections are shown in Fig. Q] (right) along 
with the relevant SED template tracks. The gap between the reddest and bluest 
K-corrections becomes exceptionally large towards high redshift (~8 mag by 
z ~ 4), which means galaxies with these different SED types are observable over 
vastly different volumes at fixed luminosity. This is a large potential source of 
selection bias if it is not accounted for in the subsequent analysis, and is one of 
the major drawbacks of optical wavelength studies at high redshift. 



Using our K-corrections, the Coe et al. ( 20061 ) photometric redshifts and a 
cosmological model with CIm = 0.3, Q\ = 0.7 and Hq = 100 km s _1 Mpc -1 we 
were able to compute absolute -B^GC-band magnitudes and sizes. The iUDF- 
BRIGHT absolute magnitude - size distribution is shown in Fig. Q] (left). Broken, 
grey lines indicate the selection limits on bright and faint galaxies for our spec- 
tral types with the reddest and bluest K-corrections. Rectangular boxes mark 
the bias free magnitude ranges of five volume limited samples selected over in- 
tervals of 0.25 in redshift from z = 0.25 to 1.5. The broken line indicates the 
Bmgc < — 18 mag cut later implemented for exposing the colour-log(n) relation. 



3. Selection Limits 

Artificial galaxy simulations are used to test our detection completeness and 
measurement reliability for galaxies in the ACS z-band image as a function of lu- 
minosit y and size. These simulations are described in detail in lCameron &: Driver 
(2007a|). In brief, we divide the luminosity-size plane into a grid and at each grid 
point generate 100 artificial galaxies using the IRAF artdata package. These 
are inserted into the i-band image and passed through our detection and mea- 
surement pipeline. The output detections, fluxes and sizes are compared to the 
input catalogue. The results of these simulations for both exponential (diffuse, 
n = 1, left-hand panels) and de Vaucouleurs (compact, n = 4, right-hand pan- 
els) profile shapes are shown in Fig. [2j Both types are readily detectable over 
the entire parameter space spanned by the observed data. However, one cannot 



Figure 1. (Left) The iUDF-BRIGHT galaxy absolute magnitude - redshift 
distribution to z = 2.5. (Right) K-corrections from observed ACS i-band to 
rest-frame MGC B-band for iUDF-BRIGHT galaxies to z = 6. 



recover accurate fluxes and half light radii beyond mean surface brightnesses of 
~27.0 mag arcsec -2 for the exponentials and only ~ 26.0 mag arcsec -2 for the 
de Vaucouleurs galaxies. The error vector diagrams indicate these mis-measured 
galaxies are erroneously recovered as faint, compact systems! Although the de 
Vaucouleurs type galaxies are easier to detect (even at faint, mean surface bright- 
nesses) due to their bright central cores, these profile shapes have extended low 
surface brightness wings, which contain a significant fraction of their total flux. 
Hence, it is more difficult to accurately recover their true magnitudes and sizes. 
This illustrates the importance of considering reliability and profile shape when 
computing selection boundaries. 



4. The Colour-Log(n) Bimodality 



Using the broadband magnitudes and Sersic indices from Coe et al. ( 20061 ) we 



are able to construct colour-log(n) diagrams for galaxies in each of our iUDF- 
BRIGHT volume-limited samples, as shown in Fig. We use our SED fits and 
K-correction software to correct from the observed filter combination straddlin 



the 40 00A break at each redshift to rest-frame Sloan (u-r) colour. iDriver et al 



(2006) have demonstrated that the local galaxy population shows a clear bi- 



modality in the (it-r) res t-log(n) plane for the Millennium Galaxy Catalogue sam- 
ple. This bimodality is most clearly seen for galaxies with abosolute magnitude 
greater than ~ — 18.0 mag in the -E^/cc-band. If fainter galaxies are included, 
these are predominately found to be blue, diffuse systems, and the blue peak 
overwhelms the distribution. Likewise, if only very bright galaxies are included, 
the bimodality is dominated by compact, red systems. Even with this cut, our 
UDF volume- limited samples do not show the colour-log(n) bimodality as clearly 
as the MGC data. This is probably due in part to small number statistics and 
in part to the large errors in UDF magnitudes, colours and photometric red- 
shifts. In addition, there appears to be a certain degree of colour evolution 
towards bluer (u-r) Tes t values at higher redshift, particularly for the more com- 
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Figure 2. (Top) Completeness, (Middle) reliability and (Bottom) error vec- 
tor diagrams for objects in the ACS UDF i-band image computed from our 
artificial galaxy simulations. The left-hand column shows the results for an ex- 
ponential (disk- like, n — 1) profile and the right-hand column shows those for 
a de Vaucouleurs (elliptical-like, n = 4) profile. The observed iUDF-BRIGHT 
galaxy distribution is marked by the black dots. 
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Figure 3. Colour-log(rt) diagrams for the MGC 0.0 < z < 0.25 sample and 
our five iUDF-BRIGHT volume limited samples from z — 0.25 to z = 1.5. 



pact galaxies. The dashed line indicates our cut at n = 1.6 used to separate 
our sample into compact and diffuse galaxy types. The n = 1.6 cut was chosen 
to lie at the saddle point between the two MGC peaks. We use a non-colour 
dependent divide as we expect a galaxy's Sersic index (i.e., its global structure) 
generally evolves on a longer timescale than its colour (i.e., its stellar population 
and recent star- formation history). 



5. Luminosity-Size Diagrams 

In Fig. |4]we present luminosity-size diagrams for galaxies in our volume limited 
samples with Bmgc < — 18 mag and 0.25 < z < 1.5, subdivided into diffuse 
(n < 1.6) and compact (n > 1.6) structural types. Regions of parameter space 
excluded by our selection limits are shaded in grey. These account for detection 
completeness, measurement reliability and the K-corrections for red and blue 
galaxies. They were constructed in the M-R e plane by com bining our artificial 
galaxy simulation results with the method of iDriven ( 19991 ). The distribution 



of galaxies occupying the remaining white-space should be free of these biases. 
The MGC local galaxy relation, also subdivided by structural type, is shown 
in the top panels and subsequently overlaid as contours. In the diffuse galaxy 
samples b, c and d where the distribution is free of bias, there is evidence of evo- 
lution towards higher surface brightnesses at high redshift and perhaps a change 
of slope in the luminosity-size relation. The compact galaxy samples are more 
tightly bound by the selection limits, but appear to show minimal evolution over 
this redshift range. 



6. Discussion 



The trends seen in our luminosity-size and colour-log(n) diagrams suggest that 
bright, diffuse (i.e., disc) galaxies are still evolving structurally from high to 
low redshift. The sense of this evolution is that the brightest disc galaxies are 
becoming larger towards low redshift. This is consistent with the picture from 
inside-out models of galaxy formation such as Truiillo and Pohlenl (2005) who 
find ~25% growth of disk galaxies since z ~ 1. Bright, compact (i.e., ellipti- 
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Figure 4. Luminosity-size diagrams for galaxies with Bmgc < — 18 mag 
and 0.25 < z < 1.5 in iUDF-BRIGHT, subdivided into diffuse (n < 1.6) and 
compact (n > 1.6) structural types. 



cal galaxies), on the other hand, seem to have finished their structural evolu- 
tion at high redshift, but continue to evolve in colour to low redshifts. This 
picture is broadly consistent with the hierarchical clustering scenario in which 
the brightest elliptical galaxies are assembled at intermediate-high redshift and 
their stellar populations passively evolve thereafter. However, it is also compat- 
ible with a scenario in whi ch spheroids form early via collapse and disc grow 
later (see lDriver et al.ll2006l ). The evolution of the disc galaxy M-R e slope sug- 
gests we are seeing more than simple passive fading and secular evolution of 
the disc galaxy population . For a more detailed investigation of these issues see 
Cameron fc Driver! (|2007bh 
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